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Exploratory  study  of  MM-wave  patch  antennas  for  strain 
measurement  and  crack  detection-Final  Report 

PI:  Haiying  Huang,  Dept,  of  Mechanical  and  Aerospace  Engineering,  University  of  Texas  at 
Arlington 

Co-PIs:  Ronald  L  Carter  and  JC  Chiao,  Dept,  of  Electrical  Engineering,  University  of  Texas  at 
Arlington 

Participating  students:  Uday  Tata  and  Shanchali  Deb 

Abstract:  The  goal  of  this  project  is  to  prove  the  strain  sensing  capability  of  a  patch  antenna. 
With  integrated  sensing  and  data  transmitting  capabilities,  the  patch  antennas  have  great 
potentials  to  serve  as  passive  wireless  strain  sensor.  The  linear  relationship  between  the  antenna 
resonant  frequency  and  the  applied  strain  was  derived  analytically.  The  procedures  for  designing 
the  patch  antenna  and  fabricating  them  using  conventional  photolithography  techniques  were 
established.  The  application  of  the  patch  antenna  for  strain  measurement  was  evaluated  by 
bonding  the  patch  antenna  to  an  aluminum  cantilever  beam  and  applying  loads  at  the  free  end  of 
the  cantilever  beam.  The  strain  sensitivity  of  the  antenna  was  estimated  to  be  15  kHz/microstrain. 
The  experimental  measurements  matched  very  well  with  the  analytical  prediction.  Noncontact 
interrogation  of  the  antenna  resonant  frequency  using  a  monostatic  radar  system  was  also 
investigated.  Research  issues  related  to  this  noncontact  interrogation  technique  is  discussed.  Two 
graduate  students,  Uday  Tata  and  Shanchali  Deb,  worked  in  this  project,  under  the  supervision  of 
the  Pis.  The  project  resulted  in  a  journal  submission,  a  provisional  patent  submission,  a  Master 
thesis,  and  a  conference  paper. 

1.  Introduction 

Structural  Health  Monitoring  (SHM)  is  a  technology  critical  for  the  safety  assurance  of 
mechanical,  aerospace,  and  civil  structures.  By  performing  autonomous  inspection  of  structures 
using  embedded  sensors,  a  SHM  system  can  provide  advanced  warnings  that  prevent  catastrophic 
structural  failures.  Strain  is  one  of  the  most  important  mechanical  parameters  acquired  by  SHM 
systems.  By  monitoring  strain  changes  in  load-bearing  structures,  the  factors  that  could  potential 
cause  structural  failures  such  as  excessive  loading,  vibration,  crack  development,  and 
environmental  aging,  etc.  can  be  detected. 

Various  strain  sensors  have  been  developed  in  the  past.  A  detailed  review  of  different  strain 
sensing  mechanisms  can  be  found  in  the  literature.  The  most  commonly  used  strain  sensor  is  the 
piezoresistive  thin  film  strain  gauge.  However,  since  SHM  systems  usually  require  a  network  of 
strain  sensors  in  order  to  obtain  a  distributed  assessment  of  the  structural  condition,  thin  film 
strain  gauges  are  not  favorable  for  SHM  systems  because  of  the  complexity  and  weight  penalty 
associated  with  the  wiring  of  these  strain  gauges.  Optical  fiber-based  strain  sensors,  on  the  other 
hand,  are  an  attractive  choice  of  sensor  for  SHM  systems  due  to  their  compact  size,  light  weight, 
remote  accessibility,  multiplexing  capability,  and  immunity  to  electromagnetic  interferences. 
Even  though  many  optical  fiber  strain  sensors  have  been  implemented  for  various  SHM 
applications,  the  inherent  shortcomings  of  optical  fiber  sensors,  e.g.  expensive  interrogation 
equipment,  limited  strain  range,  and  fragility,  etc.,  have  not  been  resolved. 

Recently,  wireless  strain  sensors  are  becoming  more  and  more  popular  for  SHM  systems.  These 
sensors,  equipped  with  embedded  microprocessors  and  wireless  communication  capabilities,  have 
the  potential  to  reduce  the  installation  and  maintenance  costs  of  distributed  sensor  networks. 
Many  of  these  wireless  sensors  require  batteries  to  be  operational.  Therefore,  their  life  spans  are 
limited  by  the  battery  life.  Passive  wireless  sensors  are  attractive  alternatives  because  they  do  not 
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require  any  power  supply.  These  passive  wireless  sensors  usually  consist  of  a  sensing  unit  and  a 
data  transmitting  unit.  The  electromagnetic  (EM)  impedance  of  the  sensing  unit  varies  with  the 
physical  measurand.  When  it  is  connected  to  an  antenna,  the  EM  impedance  change  of  the 
sensing  unit  shifts  the  resonant  frequency  of  antenna.  This  frequency  shift  is  then  remotely 
interrogated  using  inductive  coupling.  Because  the  sensing  unit  is  not  integrated  with  the  antenna, 
these  passive  wireless  sensors  are  usually  quite  large  in  size. 


In  this  project,  we  demonstrated  that  a  patch  antenna  itself  can  serve  as  the  sensing  unit  for  strain 
measurement.  With  integrated  sensing  and  data  transmitting  capability,  the  size  of  the  antenna 
sensor  can  be  reduced  significantly.  The  relationship  between  the  applied  strain  and  the  resonant 
frequency  of  the  patch  antenna  is  first  derived.  The  design,  fabrication,  and  evaluation  of  a  dual 
frequency  patch  antenna  for  strain  measurement  are  then  discussed.  The  agreement  between  the 
simulated  and  measured  strain  sensitivities  proved  the  feasibility  of  applying  patch  antennas  for 
strain  sensing. 


2.  Principle  of  operation 

The  patch  antenna,  as  shown  in  Figure  1,  consists 
of  a  layer  of  dielectric  substrate,  a  rectangular 
patch  printed  on  one  side  of  the  substrate,  and  a 
ground  plane  coated  on  the  other  side  of  the 
substrate.  The  rectangular  patch  and  the  ground 
plane,  both  made  of  conductive  metals,  form  an 
EM  resonant  cavity  that  radiates  at  a  specific 
resonant  frequency.  Based  on  the  transmission  line 
model,  the  resonant  frequency  of  a  rectangular 
patch  antenna  is  calculated  as 


Figure  1:  Patch  antenna. 
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where  c  is  the  velocity  of  light.  The  electrical  length  Le  of  the  antenna  is  defined  as  the  dimension 
of  the  metallic  patch  along  the  direction  of  the  radiation  mode.  The  effective  dielectric  constant 
sre  is  related  to  the  dielectric  constant  of  the  substrate  sr ,  the  substrate  thickness  /z,  and  the 
electrical  width  of  the  patch  we,  i.e. 
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The  line  extension  ALoc  is  calculated  from  the  effective  dielectric  constant  £re ,  the  substrate 
thickness  /z,  and  the  electric  width  we. 
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Assuming  the  antenna  is  subjected  to  a  tensile  strain  £L  along  its  electrical  length  direction,  the 
patch  width  and  the  substrate  thickness  will  change  due  to  Poisson’s  effect,  i.e. 

we  —  (1 —  u p£ l  h-(l-  us£L  )h0  .  (4) 

If  the  Poisson’s  ratios  of  the  metallic  patch  and  the  substrate  material,  up  and  are  the  same,  the 
ratio  we/h  remains  to  be  a  constant  as  the  tensile  strain  £L  changes.  As  a  result,  £re  in  equation  (2) 
is  independent  of  £L  and  ALoc  in  equation  (3)  is  proportional  to  the  substrate  thickness  h. 
Therefore,  the  resonant  frequency  in  equation  (1)  can  be  expressed  as 
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.  The  strain-induced  elongation, 


therefore,  will  shift  the  antenna  resonant  frequency.  At  an  unloaded  state,  the  antenna  frequency 
freso  is  calculated  from  the  antenna  length  Le0  and  substrate  thickness  h0,  i.e., 
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Under  a  strain  £L,  the  antenna  frequency  shifts  to 

Q 


fres  (£L  ) 


(l  +  £L  )+  C2hQ (l  “  USL  ) 


(6) 


(V) 


Combining  equations  (6)  and  (7),  we  can  establish  the  relationship  between  the  strain  eL  and  the 
frequency  shift  Afres,  i.e. 
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where  Afres=fres  -  fres0.  Analyzing  the  constant  C  indicates  that  the  strain  sensitivity  of  antenna  is 
strongly  governed  by  the  dielectric  constant  of  the  substrate  material. 


For  a  positive  strain  u  the  antenna  resonant  frequency  will  shift  to  a  lower  frequency  because  of 
the  elongation  of  its  electrical  length.  A  strain  applied  along  the  antenna’s  electrical  width 
direction,  however,  will  increase  the  antenna  resonant  frequency  due  to  the  Possion’s  effect. 
Therefore,  fres  is  more  sensitive  to  strains  applied  along  the  electrical  length  direction  than  those 
applied  along  the  electrical  width  direction.  If  the  antenna  is  fed  with  an  electric  field  polarized 
parallel  to  its  geometric  length  direction,  the  electrical  length  Le  is  the  same  as  the  geometric 
length  L  of  the  antenna.  Rotating  the  electric  field  polarization  by  90°,  however,  will  switch  the 
electrical  length  to  the  antenna’s  geometric  width  w.  By  feeding  the  antenna  with  two  orthogonal 
electric  fields,  an  antenna  that  is  sensitive  to  strains  along  both  directions  of  the  antenna  patch  can 
be  achieved. 


3.  Antenna  design 

In  order  to  feed  the  patch  antenna  with  two  orthogonal  electric  fields,  a  dual-frequency  patch 
antenna  was  designed.  For  a  dual  frequency  antenna,  the  radiation  mode  with  its  polarized  electric 
fields  parallel  to  the  geometric  length  direction  of  the  metallic  patch  is  denoted  as  the  TM10  mode 
while  the  radiation  mode  whose  electric  fields  are  along  the  geometric  width  direction  of  the 
metallic  patch  is  designated  as  the  TM0i  mode.  For  a  chosen  substrate  thickness  h  and  a  relative 
dielectric  constant  £n  the  dimensions  of  the  antenna  patch  are  calculated  from  the  resonant 
frequencies  corresponding  to  these  two  radiation  modes.  Denoting  the  resonant  frequency  for  the 
TM10  mode  as  f10,  the  geometric  length  of  the  antenna  patch  can  be  calculated  from  equation  (1) 
as 


L  = 


(9) 


Similarly,  the  geometric  width  of  the  antenna  patch  is  calculated  from  the  resonant  frequency  f0i 
of  the  TMqi  mode,  i.e. 
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where  the  line  extension  Awoc  is  calculated  from  equation  (3)  by  replacing  w  with  L.  For  a  50  pm 
thick  polyimide  substrate  with  a  relative  permittivity  of  3.4,  a  dual  frequency  patch  antenna 
operating  at  f10=  15  GHz  and  f0i  =  20  GHz  will  result  in  a  metallic  patch  with  a  length  of  5.33  mm 
and  a  width  of  4  mm. 


The  patch  antenna  can  be  fed  using  either  a  contact  or  a  non-contact  technique.  Since  this  paper  is 
focused  on  proving  the  feasibility  of  the  patch  antenna  for  strain  sensing,  the  contact  technique 
was  adopted.  In  order  to  excite  the  two  orthogonal  radiation  modes  simultaneously,  the  patch 
antenna  has  to  be  fed  at  a  proper  position.  For  the  coordinated  system  given  in  Figure  2. a,  only 
the  TM10  mode  will  be  excited  if  the  patch  antenna  is  fed  at  any  point  along  y  =  w/2.  The 
optimum  position  of  the  feeding  point  x0  is  determined  by  impedance  matching,  i.  e. , 
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where  the  calculation  of  antenna  impedance  Riny  is  given  by  Balanis  (2005).  Similarly,  only  the 
TM0i  mode  can  be  excited  if  the  antenna  is  fed  at  a  point  along  x  =  L/2.  The  feeding  position  y0 
for  TM0i  excitation  is  again  calculated  from  the  antenna  impedance  Rinx ,  i.e. 
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Feeding  the  patch  antenna  at  (x(h  y0)  excites  both  radiation  modes  and  matches  the  impedances  of 
the  patch  antenna  at  the  two  resonant  frequencies,  f10  and  f0i  (Chen  and  Wong  1996).  In  order  to 
facilitate  mechanical  testing  of  the  patch  antenna,  the  patch  antenna  was  designed  to  be  fed  using 
a  microstrip  line,  as  shown  in  Figure  2.b.  The  microstrip  line  joins  the  metallic  patch  at  (0,  y0), 
which  means  only  the  impedance  of  the  patch  antenna  for  the  TM0i  mode  is  matched.  The  length 
and  the  width  of  the  microstrip  line  were  tuned  to  achieve  the  optimum  return  losses  at  the  two 
resonant  frequencies.  The  final  design  of  the  dual  frequency  patch  antenna  is  shown  in  Figure  2.b. 
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4.  Antenna  fabrication 

The  patch  antennas  were  fabricated  on  a  flexible  polyimide  film  (Kapton  HN)  using 
semiconductor  fabrication  processes,  as  illustrated  in  Figure  3.  First,  a  50  pm  thick  Kapton  film 
was  cut  into  4”  x  4”  square  strips.  They  were  cleaned  in  ultrasonic  acetone  bath  for  a  few 
minutes.  Following  cleaning,  a  layer  of  hexamethyldisiliazane  (HMDS)  was  spin-coated  on  the 
Kapton  strips  at  a  speed  of  2000  rpm  for  30  seconds.  Positive  photoresist  SI 8 13  was  then  spin- 
coated  using  the  same  speed  and  time,  producing  a  film  thickness  of  1.5  pm  (Figure  3. a).  After 
baking  the  Kapton  strips  on  a  hotplate  at  100°C  for  one  minute,  the  antenna  pattern  was 
transferred  from  a  dark  field  mask  to  the  photoresist  by  exposing  the  photoresist  with  a  UV  light 
of  25  mW/cm2  in  intensity  for  15  seconds.  A  photoresist  developer  MF319  was  used  to  develop 
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Figure  3:  Fabrication  steps  of  the  patch  antenna. 


the  antenna  patch  (Figure  3.b).  The  developing  time  was  set  to  be  30  seconds.  After  developing, 
the  Kapton  strips  are  rinsed  with  DI  water  and  dried  using  dry  nitrogen  air.  A  copper  film  of  1  pm 
in  thickness  was  then  deposited  on  the  Kapton  strips  using  thermal  evaporation  (Figure  3.c).  The 
final  fabrication  involved  soaking  the  Kapton  strips  in  pure  acetone  for  3  hours  to  remove  the 
developed  photoresist  and  lift  off  the  metal  on  top  (Figure  3.d).  Since  we  will  bond  the  antenna  to 
a  metallic  structure,  no  ground  plane  was  fabricated  on  the  substrate.  The  fabricated  patch 
antenna  is  shown  in  Figure  4. 


5.  Mechanical  testing  of  the  patch  antenna 

Mechanical  testing  of  the  patch  antenna  for  strain  measurement  was  carried  out  using  a  cantilever 
beam,  as  shown  in  Figure  5. a.  The  antenna  sensor  was 
bonded  on  an  aluminum  cantilever  beam  using  a 
conventional  strain  gauge  adhesive.  The  microwave 
signal  supplied  by  a  network  analyzer  was  fed  to  the 
antenna  sensor  through  an  SMA  connector.  The  signal 
reflected  by  the  antenna  sensor  was  routed  back  to  the 
receiver  of  the  network  analyzer,  from  which  the  antenna 
radiation  parameters,  i.e.  the  Su  parameter,  can  be 
determined.  Weight  was  added  at  the  free  end  of  the 
cantilever  beam  at  an  increment  of  5  lb.  Figure  4:  Fabricated  patch  antenna. 
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Figure  5:  Experimental  set  up  for  strain  measurement;  (a)  loading  of  the  cantilever  beam;  (b) 
width-direction  elongation;  (c)  length-direction  elongation. 


In  order  to  stretch  the  antenna  along  its  width  direction,  the  patch  was  oriented  with  its  microstrip 
line  perpendicular  to  the  length  direction  of  the  cantilever  beam  (see  Figure  5.b).  An  SMA 
connector  was  mounted  on  the  edge  of  the  cantilever  with  its  pin  touching  the  end  of  the 
microstrip  line.  A  solder  ball  was  placed  at  the  end  of  the  SMA  tip  to  ensure  a  sound  electrical 
connection  between  the  SMA  tip  and  the  antenna  microstrip  line.  In  order  to  stretch  the  antenna 
along  its  length  direction,  however,  the  antenna  has  to  be  oriented  with  its  microstrip  line  parallel 
to  the  length  direction  of  the  cantilever  beam  (see  Figure  5.c).  The  only  feasible  way  to  feed  the 
patch  antenna  was  then  to  vertically  hold  the  SMA  connector  and  press  the  tip  of  the  SMA 
connector  onto  the  end  of  the  antenna  microstrip  line. 

6.  Results  and  discussions 

The  design  of  the  dual  frequency  antenna  was  verified  using  a  commercial  EM  simulation  tool, 
Sonnet  Pro  vll.5.  The  simulated  Sn  curve  of  the  dual  frequency  antenna  is  shown  in  Figure  6. 
The  two  resonant  frequencies  were  found  to 
b q  f jo  =  14.9  GHz  and  f0i  =  19.6  GHz.  The 
return  loss  at  the  f0i  is  higher  than  that  at  f10 
because  of  better  impedance  matching  for  the 
TM0i  mode.  The  Sn  parameters  of  the  patch 
antenna  under  zero  loading  were  measured 
and  compared  with  the  simulated  results.  The 
measured  patch  antenna  displayed  two 
resonant  frequencies  in  its  Sn  curve.  The 
length-direction  resonant  frequency  f10  was 
measured  at  17.2  GHz  while  the  width- 
direction  resonant  frequency  f0i  was  measured 
at  20.5  GHz.  The  2  GHz  discrepancy  between 
the  designed  frequency  and  the  measured 
frequency  for  the  TM10  mode  was  probably 
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Figure  6:  Comparison  of  simulated  and 
measured  antenna  SI  1  curves. 
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Figure  7:  Shifts  of  Sn  curves  under  width-direction  strains;  (a)  overview;  (b)  shifts  of f10 
frequency-enlarged  view;  (c)  shifts  of  the  f0i  frequency-enlarged  view. 


due  to  the  additional  capacitance  introduced  by  the  SMA  connector. 


The  patch  antenna  was  tested  by  applying  strains  along  its  geometric  width  direction  first.  The 
load  was  increased  at  an  increment  of  5  lb  until  a  total  load  of  30  lb.  At  the  antenna  position,  5  lb 
of  load  will  lead  to  0.084%  of  strain.  The  network  analyzer  was  set  to  acquire  1600  points  from 
16  GHz  to  21.5  GHz.  An  overview  of  the  measured  Sn  curves  under  different  strain  levels  is 
shown  in  Figure  7. a.  Zooming  in  the  regions  near  the  foi  and  f10  frequencies  (see  Figures  7.b  and 
7.c)  demonstrated  the  parallel  shifts  of  these  two  frequencies.  The  f10  frequencies  shifted  to  the 
right  as  the  strains  increased  while  the  f0i  frequencies  shifted  to  the  left.  This  is  consistent  with  the 
analytical  prediction.  The  percentages  of  frequency  shifts  for  these  two  resonant  frequencys,  both 
calculated  and  measured,  were  plotted  with  respect  to  the  applied  strains,  as  shown  in  Figure  8. 
Excellent  linearity  was  achieved  for  both  the  calculated  and  measured  data.  The  slopes  of  the  fit 
lines  represent  the  strain  sensitivity  of  the  respective  frequency.  For  the  f0i  frequency,  the 
analytical  model  predicted  a  slope  of  0.98,  which  means  a  1%  strain  would  shift  the  f0i  frequency 
by  0.98%.  The  slope  of  the  measured  data  is  0.81,  slightly  lower  than  that  of  the  calculated  data. 
This  is  probably  due  to  the  shear  lag  introduced  by  the  bonding  layer  and  the  Kapton  substrate. 
Expressing  the  frequency  shifts  in  absolute  value,  the  strain  sensitivity  of  the  f10  frequency  was 
estimated  to  be  16.4  kHz/microstrain.  Assuming  the  network  analyzer  can  achieve  a  10  Hz 
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Strain  along  width  direction  (%) 


Figure  8:  Comparison  of  analytical  and  measured  strain  sensitivities; 
(a)  fio  frequency;  (b)  f0i  frequency. 


frequency  resolution,  the  strain  resolution  of  the  antenna  sensor  is  estimated  to  be  6.1e-4 
microstrains.  The  sensitivity  of  the  f0i  frequency  to  the  width  direction  strains  is  much  smaller. 
With  a  slope  of  0.18,  the  strain  sensitivity  in  absolute  value  is  3.3  kHz/microstrain,  one  fifth  of 
the  strain  sensitivity  of  the  f10  frequency. 

When  the  antenna  was  strained  along  its  length  direction,  only  the  f10  frequency  can  be  observed 
at  the  network  analyzer  due  to  vertical  feeding.  Therefore,  the  network  analyzer  was  set  to  acquire 
400  data  points  from  17.4  GHz  to  18.1  GHz.  Because  the  bending  of  the  cantilever  beam  causes 
the  antenna  to  slide  away  from  the  SMA  connector,  the  antenna  had  to  be  mounted  close  to  the 
fixed  end  of  the  cantilever  beam.  Applying  a  load  of  5  lb  at  the  free  end  of  the  cantilever  beam 
will  produce  a  strain  of  0.1%  at  the  antenna  position.  The  parallel  shifts  of  the  Sn  curve  under 
various  strains  are  shown  in  Figure  9. a.  Again,  we  achieved  an  excellent  linearity  between  the 
shifts  of  the  f10  frequency  and  the  length-direction  strains  (see  Figure  9.b).  Moreover,  the  slope  of 
the  measured  data  matched  very  well  with  the  theoretical  predictions.  The  sensitivity  of  the  f10 
frequency  with  respect  to  the  length-direction  strain,  around  17.2  kHz/microstrain,  is  slightly 
higher  than  that  of  the  foi  frequency  when  the  antenna  was  under  width  direction  elongation. 


17.4 


17.6 


17.8 


Frequency  (GHz) 
18  18.2 


Strains  along  length  direction  (%) 
0.1  0.2  0.3  0.4  0.5 


Figure  9:  (a)  Shifts  of f10  frequency  under  length  direction  strains;  (b)  comparison  of 
analytical  and  measured  strain  sensitivities. 
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7.  Discussion  of  non-contact  measurement  of  antenna  resonant  frequency 

The  antenna  resonant  frequency  can  be  measured 
via  a  non-contact  reader,  based  on  the  principle  of 
backscattering,  as  shown  in  Figure  10.  The  reader 
serves  as  both  the  transmitter  and  the  receiver.  As  a 
transmitter,  the  reader  sends  an  incident  EM  wave 
toward  the  patch  antenna.  Upon  intercepted  by  the 
patch  antenna,  the  incident  EM  wave  is  scattered 
back  by  the  patch  antenna  if  the  frequency  of  the 
incident  wave  matches  the  resonant  frequency  of 
the  patch  antenna.  The  backscattered  EM  wave  is 
then  received  by  the  reader.  By  performing 
frequency  analysis  of  the  backscattered  EM  signals, 
the  antenna  radiation  parameters  (resonant 
frequency,  return  loss,  etc.)  can  be  measured. 


Figure  10:  Non-contact  interrogation  of 
antenna  resonant  frequency 


Assuming  a  horn  antenna  is  employed  as  the  antenna  reader  and  it  is  placed  at  a  distance  R  from 
the  antenna  sensor,  the  power  of  the  backscattered  signal  received  by  the  antenna  reader  can  be 
calculated  from  the  Friis  Transmission  Equation,  i.e. 


p,gIg;a4d 

r  (4  xR)4 


(12) 


where  P,  is  the  transmitting  power,  Gh  is  the  gain  of  the  horn  antennas,  Gs  is  the  gain  of  the 
antenna  sensor,  X  is  the  wavelength,  and  D  is  the  insertion  loss  of  the  sensor.  The  minimum 
detectable  power  at  the  receiver  is  calculated  as 

Rmin  =  SNR  xkTxBx  NF  ,  (13) 


where  SNR  is  the  signal-to-noise  ratio,  kT  is  the  thermal  energy,  B  is  the  bandwidth  of  the 
receiver,  and  NF  is  the  noise  Figure  of  the  receiver.  The  maximum  distance  at  which  the  antenna 
sensor  can  be  detected  is  therefore 


^max  , 

4ft 


DPtG2hGl 


\  1  /  4 


SNR  xkTxBx  NF 


For  a  transmitting  power  of  Pt=  1  mW,  Gh 
=  316  (25  dBi  gain),  Gs  =3.98  (6  dBi  gain), 
X  =  15  mm  (RE  =  20  GHz),  SNR  =  100,  kT 
=  4xl0"21  J,  B  =  10  kHz,  NF  =  1,  and  D  = 
0.1,  the  estimated  maximum  detectable 
distance  is  16.8  meters. 


(14) 


A  waveguide  was  originally  used  as  the 
antenna  reader.  As  shown  in  Figure  1 1,  the 
waveguide  is  mounted  on  a  three-axis 
translation  stage  for  proper  alignment 
between  its  opening  and  the  patch  antenna 
below  it.  A  MATLAB  program  was 
developed  to  time-gate  the  time  domain 
antenna  Sll  signal.  By  convert  the  time- 


Figure  11:  Remote  interrogation  of  antenna 
resonant  frequency  using  a  waveguide 
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Figure  12:  Time  gated  Sll  signal;  (1)  time-domain;  (2)  frequency  spectrum 


gated  Sll  signal  back  into  the  frequency  domain,  the  antenna  resonant  frequency  can  be 
determined.  Figure  12. a  shows  a  typical  gated  Sll  signal  in  time  domain.  The  spectrum  of  the 
time-gated  Sll  signal  matched  well  with  the  spectrum  of  the  SI 2  signal.  However,  the  shifts  of 
the  Sll  signal  were  not  consistent  with  the  applied  strain  when  it  was  tested  using  the  cantilever 
beam.  The  challenge  facing  in  using  the  waveguide  as  the  antenna  reader  is  that  the  spectrum  of 
the  Sll  signal  is  highly  sensitive  to  the  waveguide-antenna  alignment.  A  small  movement  in  the 
waveguide  position  will  disturb  its  alignment  with  the  patch.  In  addition,  the  waveguide  needs  to 
be  placed  very  close  to  the  patch  since  it  has  a  low  directivity.  We  have  also  tried  to  use  a  horn 
antenna  as  the  antenna  reader.  However,  the  horn  antenna  usually  has  a  large  far-fleld  distance, 
typically  2.5  meters  at  20  GHz.  When  the  antenna  is  placed  at  the  far-field,  the  beam  size  of  the 
horn  antenna  is  very  big.  Even  with  time  gating,  the  backscattering  of  the  patch  antenna  is  of  low 
signal  strength  compared  to  the  strength  of  background  cluttering.  Because  of  the  small  signal-to- 
noise  ration,  it  is  difficult  to  determine  the 
frequency  of  the  backscattered  signal. 

It  was  decided  that  a  bi-static  radar  system  will 
help  alleviating  this  problem.  As  shown  in 
Figure  13,  employing  two  horn  antennas,  one  for 
transmitting  and  one  for  receiving,  gives  more 
flexibility  in  the  alignment.  By  properly  align  the 
polarization  of  the  two  horn  antennas,  the 
background  cluttering  can  be  reduced 
significantly.  Due  to  the  long  lead  time  for  the 
second  horn  antenna,  the  testing  of  the  bi-static 
radar  system  for  antenna  resonant  frequency 
interrogation  will  be  carried  out  in  the  future. 

8.  Conclusions 

This  report  presented  the  analysis,  design,  fabrication,  and  validation  of  a  patch  antenna  for  strain 
measurement.  Both  the  theoretical  analysis  and  the  experiment  results  demonstrated  that  the 
resonant  frequency  of  a  patch  antenna  is  sensitive  to  strains  applied  along  the  direction  of  its 
corresponding  radiation  mode.  Since  the  antenna  resonant  frequency  can  be  remotely  interrogated 
using  field  coupling  or  antenna  backscattering,  the  patch  antennas  have  great  potentials  to  serve 
as  passive  wireless  strain  sensor.  The  development  of  a  wireless  radar  system  that  can  interrogate 
the  antenna  resonant  frequency  remotely  is  under  way. 


/ 


Figure  13:  Bi-static  radar  system  for 
antenna  resonant  frequency  interrogation 
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Appendix  A:  List  of  publications 

1.  U.  Tata,  H.  Huang,  R.  Carter,  JC  Chiao,  “Exploiting  patch  antenna  for  strain  measurement”, 
submitted  to  Measurement  Science  and  Technology 

2.  U.  Tata,  H.  Huang,  R.  Carter,  JC  Chiao,  “Passive  Wireless  Sensor  for  Strain,  Temperature, 
Crack,  and  Fatigue  Measurement”,  US  provisional  patent  application 

3.  U.  Tata,  “Study  of  patch  antennas  for  strain  measurement”,  Master  thesis,  University  of 
Texas  at  Arlington,  Aug.  2008  (available  from  UTA  library) 

4.  U.  Tata,  H.  Huang,  R.  Carter,  JC  Chiao,  “Bio-inspired  sensor  skins  for  structural  health 
monitoring”,  to  be  presented  at  the  ASME  Conference  on  Smart  Materials,  Adaptive 
Structures  and  Intelligent  Systems,  October  28-30,  2008,  Ellicott  City,  Maryland,  USA 

Appendix  B:  Executive  summary  of  Tata’s  Master  thesis 

An  innovative  method  of  measuring  strain  using  a  patch  antenna  is  investigated  in  this  thesis  to 
overcome  the  limitations  of  existing  wireless  strain  sensing  technologies.  The  patch  antenna  is 
made  of  a  thin  sheet  of  low-loss  insulating  material,  called  the  dielectric  substrate.  The  antenna 
pattern,  i.e.  a  metallic  patch,  is  printed  on  one  side  of  the  substrate.  A  ground  plane  is  coated  on 
the  opposite  side  of  the  dielectric  substrate.  The  metallic  patch  and  the  ground  plane  form  an 
electro-magnetic  (EM)  cavity.  When  fed  with  electromagnetic  signal,  the  patch  antenna  radiates 
at  a  resonant  frequency.  The  antenna  resonant  frequency  depends  on  its  electrical  dimensions. 
Strain  changes  the  dimensions  of  patch  antenna,  resulting  in  a  shift  in  the  resonant  frequency  of 
patch  antenna.  First,  single  frequency  and  dual  frequency  antennas  are  designed  using  the 
transmission  line  model.  The  antenna  design  is  then  confirmed  using  an  Electromagnetic  (EM) 
simulation  tool,  Sonnet  11.5.  Because  the  dual  frequency  antenna  has  two  fundamental 
frequencies  f0io  and  f0io  corresponding  to  the  antenna’s  geometrical  length  and  width  respectively, 
it  is  sensitive  to  strains  applied  along  both  directions.  The  single  frequency  antenna,  however,  has 
only  one  fundamental  frequency  f0w .  Therefore,  it  is  only  sensitive  to  strains  applied  along  the 
antenna  length  direction.  The  effect  of  strain  along  the  width  and  length  direction  of  the  dual 
frequency  antenna  and  the  width  direction  of  the  single  frequency  antenna  on  the  antenna 
resonant  frequency  were  simulated  using  Sonnet  1 1.5  and  the  strain  sensitivity  of  the  antenna  was 
calculated  for  each  case.  The  antennas  are  fabricated  on  a  flexible  Kapton  substrate  using 
conventional  micromachining  techniques.  The  effect  of  strain  on  antenna  resonant  frequency  was 
verified  experimentally  by  bonding  the  antenna  onto  a  cantilever  beam  and  applying  load  at  one 
end  of  the  cantilever  beam.  The  experimental  and  simulated  results  are  in  good  agreement  with 
each  other. 

Appendix  C:  List  of  software  program 

1.  MATLAB  program  to  simulate  the  effect  of  strain  on  the  resonant  frequency  of  a  patch 
antenna  based  on  transmission  line  model 


%  Calculate  the  frequency  of  microstrip  patch  antenna 
%  from  Bhartia,  Rao  &  Tomar,  "Millimeter-wave  microstrip  and  printed 
%  circuit  antennas ", Artech  House,  Boston,  1991,  P10-11 
%  First  created:  10/28/2006 
clear;  clc; 

w0=0.004;  %  width  of  patch,  in  meter 
L0=0 . 00533  %  length  of  patch  in  meter 
t0=50e-6;  %  thickness  of  substrate,  in  meter 
Epsilon=3.4;  %  dielectric  constant  of  the  substrate 

strain= [ 0 : 0 . 01 : 1 ] *le-2 ;  %  strain  of  antenna 

mu_c=0.3;  %  poissons  ratio  of  the  conductor 
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mu_s=0.3;  %  poisson's  ratio  of  the  substrate 

L= ( 1+strain) *L0 ;  %  length  of  patch  after  strain 

w=  ( l-mu_c*strain)  *w0 ;  %  width  of  patch  after  strain 

t= ( l-mu_s*strain) *t0 ;  %  thickness  of  the  substrate  after  strain 

Rwt=w . /t ; 

%alpha=l+l/49*log ( (RwtA4+ (l/52*Rwt) A2) / (RwtA4+0 . 432) ) +1/18 . 7*log (1+ (1/18 
. 1  *Rwt ) . A3) ; 

Eref f= (Epsilon+1) /2+ (Epsilon-1) /2* (1  +  10. /Rwt)  . A  (-0 . 5) ;  %  effective 

dielectric  constants  of  substrate 

dL=0 .412*  (Eref f +0.3)  .  * (Rwt+0.264)  ./  (Eref f-0 . 258 )  ./ (Rwt+0.813)  .  *t;  % 

line  extension 

f r=2 . 998e8/2 . /sqrt (Eref f ) . / (L+2*dL) ;  %  resonant  frequency  of  the 

antenna 

plot ( strain*100 ,  fr/le9);  title  (’ Frequency  vs.  strain '); xlabel (' strain 
(%) 1 ) ; ylabel ( 1 R. F  (GHz ) 1 ) ;  grid  on; 

(max (fr) -min (fr) ) / (max (strain) *le6) 

2.  MATLAB  program  to  evaluate  the  shift  of  the  SI  1  curve 

%interpolation  of  the  frequency  values  for  a  given  set  of  returnloss 
values 

%  R  is  resonant  frequency  values 

%  Zero,  five,  ten,  fifteen,  twenty, twentyfive  and  thirty  are  return 
%  loss  values  for  all  the  loads 
%  xi  chosen  set  of  returnloss  values 

y=  [  ]  ; 

zero=load ( ' zero . txt 1 ) ;  %  load  the  return  loss  values  from  the  txt  file 

f ive=load ( ' five . txt 1  )  ; 

ten=load ( 1  ten . txt 1  )  ; 

f if teen=load ( 1  fifteen . txt 1  )  ; 

twenty=load ( 1  twenty . txt ' ) ; 

twentyf ive=load ( 1 twentyfive . txt 1 ) ; 

thirty=load ( ' thirty . txt 1 ) ; 

Reload ( ' R . txt 1 )  ; 
for  xi=-14 : -0 . 1 : -15 

Rxi=f_R_xi (thirty,  R,  xi) ;  %  function  to  interpolated  the  returnloss 
values 

y= [y, Rxi ] ; 

end 

y  (  :  ) 

function  Rxi=f_R_xi (x,  R,  xi) 

%  function  Rxi=R_kxi (x,  R_meas,  xi) 

%  find  the  value  R  at  xi  if  the  relationship  between  x  and  R  is  given 
%  by  R  at  location  x.  x  and  R  are  sort  in  ascend  order 
%  Inputs: 

%  x-sampled  locations 
%  R-value  at  sampled  locations 
%  xi-specified  locations 
%  Outputs : 

%  Rxi-value  at  xi  derived  from  x,  R 
N=length (xi ) ; 
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for  ii=l : N 

index=find (x==xi (ii)  )  ; 

if  (length (index) ==1 ) 

Rxi ( ii ) =R ( index)  ; 

else 

index=find (x<xi (ii)  )  ; 
p_l=index (length (index)  )  ; 

x _ 1— x (p_l) ; 

index=find (x>xi (ii)  )  ; 
p_h=index ( 1 )  ; 
x_h=x (p_h) ; 

Rxi ( ii ) =R (p_l ) + (xi ( ii ) -x_l ) / (x_h-x_l ) * (R (p_h) -R (p_l ) ) ; 

end 

end 

[M,  N]  =size  (R)  ; 

if  N==l 

Rxi=Rxi . 1 ; 

end 


3.  MATLAB  program  for  converting  the  SI  1  curve  to  time-domain  signal  time  gating  the  time- 
domain  signal  and  converting  the  time-gated  signal  back  to  frequency  domain 

function  [t,  R,  fk,  Rf ft ] =f_2port_S (fn) 

%  this  program  process  the  S  parameter  obtained 
%  by  a  two  port  measurement 

%  IFFt  is  used  to  convert  the  S  parameters  to  time  domain  signals 
%  the  data  were  loaded  from  fn.mat 
%  f :  frequency 

%  S :  S  parameter  in  complex  number 

%  Created  on  04/07/08 


C= [ 1  load  ',  fn,  '.mat']; 
eval (C) ;  %  load  the  data 

f 0= [ 0 : f ( 2 ) —  f (1)  : min ( f )  ]  .  '; 
f  =  [  f  0 ;  f  ]  ; 

S12=[zeros  (size  (fO) ) ; S 1 2 ] ; 

figure;  plot(f,  abs (S12 )); title ( 1 S  parameter');  xlabel('Freq  (GHz)'); 
ylabel  (  ' S  (dB)  ' ) ;  grid  on 

Nfft=2A15; 


%  convert  the  S  data  back  to  time  domain 
[t,  R]=f_ifft_t (S,  f (2 ) -f ( 1 ) ,  Nfft); 

figure; plot (t,  R) ;  title ('time  domain  signal  S');  grid  on 
tl=input (' input  low  cut-off  signal:  '); 
th=input (' input  high  cut-off  signal:  '); 


%  window  the  antenna 
index=f ind (t>tl) ;  % 

tc-t ( index) ; 

Rc=R ( index)  ; 
index=f ind ( tc<th) ;  % 
tc=tc ( index)  ; 


signal 

cut-off  the  signals 


high  cut-off 


at  the  begining 
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Rc=Rc ( index) ; 

%Rc=Rc . *hanning (length (Rc) ) ; 

Ntc=length (tc) ; 

[fk,  Rf ft] =f_f ft_f k (Rc,  1/ (tc  (2 ) -tc  (1) ) ,  Nfft); 

figure;  plot  (tc,  Rc) ; title (' Windowed  time  domain  signal');  xlabel('time 
(ns) '  )  ;  ylabel ( 'Amp' ) ;  grid  on; 

figure;  plot(fk,  Rf ft/max (Rf ft )) ;  title  (' S  from  windowed  time  domain 
signal');  grid  on;axis([0  50  0  1]); 
xlabel('Freq  (GHz)');  ylabel (' amp ') ; 

return; 

function  [t,  Rif ft] =f_if ft_t (Ref ft,  fs,  Nfft) 

%  this  function  convert  a  frequency  domain  signal  to  time  domain  signal 
%  Inputs : 

%  Ref ft:  frequency  domain  signal  in  complex 

%  fs:  frequency  sampling  rate 

%  Nfft:  number  of  inverse  FFT 

%  Outputs : 

%  t:  time 

%  Rif ft:  signal  in  time  domain,  obtained  from  inverse  fft 

Rifft=ifft (Rcfft,  Nfft) ;  %  digital  inverse  fft 

Rif ft=real (Rif ft (1 : floor (Nf ft/2 ))) ;  %  the  real  part  only 

t=l/fs* [0 : floor (Nfft/ 2) -1] /Nfft; 

return 

function  [fk,  Rf f t ] =f_f f t_f k (R,  fs,  Nfft) 

%  function  [fk,  Rf ft ] =f ft_f k (R,  fs,  Nfft) 

%  find  the  FFT  of  R  as  a  function  of  frequency  fk 
%  Inputs: 

%  R:  signal 

%  fs:  sampling  frequency 

%  Nfft:  #  of  fft 

%  Outputs : 

%  fk:  frequency 

%  Rfft:  real  value  of  fft (R) 

Rf f t=abs (f f t (R,  Nfft));  %  fft  of  signal  R 
Rf f t=Rf f t ( 1 : floor (Nfft/2)  )  ; 
fk=fs* [0:floor (Nfft/2) -1] /Nfft; 
return 
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